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ABSTRACT 

Context. With the legacy of Spitzer and current advances in (sub)mm astronomy, a considerable number of so-called 'transitional' 
disks has been identified which are believed to contain gaps or have developped large inner holes, some filled with dust. This may 
indicate that complex geometries may be a key feature in disk evolution that has to be understood and modelled correctly. The disk 
around LkCa 15 is such a disk, with a large gap ranging from ~5 - 46 AU, as identified by Espaillat et al. (2007) using 1 + lD radiative 
transfer modelling. To fit the spectral energy distribution (SED), they propose two possible scenarios for the inner (<5 AU) disk - 
optically thick or optically thin - and one scenario for the outer disk. 

Aims. We use the gapped disk of LkCa 15 as a case in point to illustrate the importance of 2D radiative transfer in transitional disks 
by showing how the vertical dust distribution in dust-filled inner holes determines not only the radial optical depth but also the outer 
disk geometry. 

Methods. We use MCMax, a 2D radiative transfer code with a self-consistent vertical density and temperature structure, to model the 
SED of LkCa 15. 

Results. We identify two possible geometries for the inner and outer disk that are both dilferent from those in Espaillat et al. (2007). 
An inner disk in hydrostatic equilibrium reprocesses enough starlight to fit the near infrared flux, but also casts a shadow on the inner 
rim of the outer disk. This requires the outer disk scale height to be high enough to rise out of the shadow. An optically thin inner 
disk does not cast such a shadow, and the SED can be fitted with a smaller outer disk scale height. For the dust in the inner regions 
to become optically thin however, the scale height would have to be so much higher than its hydrostatic equilibrium value that it 
effectively becomes a dust shell. It is currently unclear if a physical mechanism exists which could provide for such a configuration. 
Conclusions. We find that the radial optical depth of dust within the inner hole of LkCa 15 is controlled by its vertical distribution. If 
it turns optically thick, the outer disk scale height must be increased to raise the outer disk out of the inner disk's shadow. 

Key words. Stars: pre-main sequence - Stars: individual: LkCa 15 - planetary systems: protoplanetary disks - radiative transfer - 
circumstellar matter 

1. Introduction thin. Additional observations show that this dust is most likely 

present in the form of a small optically thic k inner disk at the 
Protoplanetary disks are the main sites of planet formation, dust evaporating radius (lEspaiUat et allEOQ^ hereafter Esp08). 
Within them, small dust grains grow from micron to millime- j^g inbetween is mostly empty 

ter sizes and settle to the midplane where they are beheved to j^is paper we will model both geometries proposed by 

eventually form planetary systems. Sufficiently massive planets ^^^q-j ^^^^^^ ^ f^jj 2d radiative transfer code and focus on the 



5^ locally deplete the disk of gas and dust and may open a gap or vertical structure of the inner and outer disk. Our modelling ap- 

5^ , enlarge the inner hole. Recently, considerable numbers of so- ^^.^^^^ differs from that of Esp07 in one important point: We 

called transitional disks have been identified on the basis of their combine the several disk components in one model before we 

low near infrared excess ('Nanta et aU2007|), indicating that their perform the radiative transfer and create the spectral energy dis- 

inner regions may be (partially) depleted of dust. tribution (SED), whereas they create the SED after perform- 

LkCa 15 is such a ti-ansitional disk, but with a more com- ing the radiative transfer on each seperate component. Although 

plex geometry: It has an inner disk, a gap and an outer disk, their approach is valid if the seperate components do not influ- 

The oute r disk, seen at an inclination of 51°, extends from 46 to ence each other, we will show that at least for LkCa 15 this is 

800 AU (Pietu et al. 2006). The Spitzer data show that t he inner not the case. 
hole i s not devoid of dust between 0.1 and 5 AU (Espaill at et al.l 
l2007l hereafter Esp07). On the basis of 1h-1D RT modelfing, 
Esp07 identify two possible scenarios for the inner regions: one 

where the dust is optically thick and one where it is optically xhe disk model we use here is MCMax (iMin et al.ll2009l) . a 2D 

Monte Carlo radiative transfer code with a selfconsistent verti- 

Send offprint requests to: G.D.Mulders, e-mail: muldersOuva.nl cal structure. The Monte Carlo radiative transfer simulation is 



2. Disk model 



2 



Mulders, Dominik, Min: Full two-dimensional radiative transfer modelling of the transitional disk LkCa 15 (RN) 



based on that of iBiorkman &Woodl (120011) and calculates the 
dust temperature for a given density structure. It also calcu- 
lates the vertical structure by directly integrating the equation 
of hydrostatic equilibrium using the Monte Carlo temperature. 
Because the vertical structure calculation requires the tempera- 
ture in the midplane to be accurately determined, MCMax treats 
those reg ions with an analytical diffusion approximation (see 
iMin et al . (2009) for details). By iterating between the Monte 
Carlo simulation and the hydrostatic equilibrium calculation a 
self-consistent solution is found. This typically takes four or five 
iterations. MCMax also includes a modified random-walk ap- 
proximation that speeds up the calculations in regions with ex- 
treme optical depth s by making mu ltiple inte raction steps in a 
single computation (iMin et al. I ( 120091) . see also lFleck & Canfieldl 
( |1984|) ). The typical run time is less than three minutes on a 
single core CPU for a model with a radial optical depth in the 
midplane of Ti^m = 10^ using six iterations with 10^ photon 
packages each. 

Because models with well-mixed dust and gas in vertical 
hydrostatic equilibrium tend to overpredict the disk's surface 
height and far infrared flux, dust settlin g and grain growth ar e 
necessary ingredients to fit SEDs (e.g. iD'Alessio et ani2006h . 
For the present study we will describe them with two parameters 
that lower the disk surface: The mass fraction of large grains 
fmm and a sedimentation parameter ^F^n. The sedimentation pa- 
rameter Tall describes the settling of small dust particles and 
is implemented by reducing the dust scale height with respect 
to the selfconsistent gas scale height: HjustCz) = *PaiiHgas(z). A 
similar global sed imentation parameter has also been used by 
iRatzka et al.l (12007 ). Both parameters have a different effect on 
the surface shape and SED and are not degenerate. The sedimen- 
tation parameter is global and reduces the scale height in both the 
inner and outer disk and therefore scales the SED accordingly. 
The fraction of large, mm-sized grains only affects the surface 
height in the outer disk: Increasing fm,n reduces the flaring angle 
of the outer disk and mostly affects the far infrared flux. 

We will use a bimodial size distribution here, consisting of 
small (O.ljum) and big (2 mm) grains. For the dust opacity of 
the small part icles we use that of .1 fim grains, consisting of 
80% silicate (Dorschner et alJ [l995h and 20% amorphous car- 
bon dPreibisch et al.. 1993) . We do not m odel them as com pact 
spheres, but as irregular shaped particles (iMin et alJl2005h . The 
opacity of the big particles is modelled using 2 mm sized irreg- 
ular shaped silicate particles. The total opacity is determined by 
the fitting parameter fmni, the fraction of mass in big particles. To 
settle the big grains to the midplane, we introduce a stratification 
parameter ^Fbig that reduces the scale height of large grains with 
respect to that of small grains: Hbig(z) = *l'bigHdu,st(z). 

For the stellar photosphere we use a main-sequence Kurucz 
model wit h the stellar mass, luminos ity and effective temper- 
ature from iKenvon & HartmannI (1199 6) (Table [TJ. This photo- 
sphere is fitted to the de-reddened photometry between 0.5 fim 
and 2 fim by scaling the distance. 

3. Results 

For each of the geometries (Fig. [TJ we follow a different ap- 
proach. In the optically thick case, we model the inner disk first 
to constrain the size of the shadow cast on the outer disk and pro- 
ceed by modelling the outer disk. Note that although we change 
the fit parameters for one component at a time, we always per- 
form the radiative transfer including all components. In the op- 
tically thin case, the inner disk produces no shadow and we can 
model the outer disk first. We then proceed by including an inner 




b) Optically thin dust shell 




Fig. 1. Sketches (not to scale) of the two possible geometries for 
LkCa 15 from this paper (top panel) and - for comparison - that 
from Esp07 (bottom panel). 

Depicted are: a) A model with a 1 AU optically thick inner disk 
and an outer disk that rises out of its shadow, b) A model with 
a 5 AU optically thin spherical dust shell and a flatter outer disk 
with a fully illuminated inner rim. c) The model proposed by 
Esp07, a combination of the inner disk from (a) with the fully 
illuminated outer disk of (b). 



disk so we can explore the effect of the vertical dust distibution 
on the radial optical depth. 

3. 1 . Optically thick inner disk 

In the optically thick case, the dust mass in the inner disk is 
hard to constrain from the SED. Therefore we use the same sur- 
face density powerlaw (p=l) as for the outer disk, which yields 
'E(IAU) = 560 g/cm^, assuming a dust-to-gas ratio of 1:100. We 
fix the inner radius at the dust evaporating radius at 0. 1 AU and 
fit an outer radius for the inner disk of 1 .0 AU, yielding a dust 
mass for the inner disk of 4 x 10 M©. Our best fit (Fig. |2]l does 
not require any sedimentation, and does not constrain stratifica- 
tion or big particle fraction, for which we will use the fit values 
from the outer disk. 

This choice of inner disk parameters is not unique, but the 
size of the shadow cast on the outer disk is. The height at the 
outer edge of the inner disk (z/r - ~ 0.1, Fig. |3]l sets the total 
amount of energy reprocessed by the inner disk, which is con- 
strained by the SED. The size of the shadow is therefore approx- 
imately equal for all solutions and independent of the exact inner 
disk geometry. 
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Fig. 2. Spectral Energy Distribution of LkCa 15 (diamonds) with 
Icr errorbars. Overplotted in grey is the stellar photosphere. The 
top panel shows our best-fit disk model with an optically thick 
inner disk (Fig.[T]a) 

The solid line in the bottom panel shows our best fit with an 
optically thin dust shell of 1 x 10"" Mq (Fig.[I]b). The dot- 
ted line shows the same outer disk combined with an inner disk 
in hydrostatic equilibrium of 2 x 10""^ M© that is optically 
thin in the vertical direction (Fig [1] c). It reproduces the near 
and mid infrared flux, but is not radially optically thin. It casts 
a shadow on the outer disk and underpredicts the far infrared 
flux by almost a factor of ten. The dashed line shows the same 
model with a very flat (z/r ~ 0.01) inner disk. This minimizes 
the size of the shadow, but such an inner disk does not reprocess 
enough of the stellar radia tion to fit the flux shortward s of 20 
urn. References - (UB VRI) iKenv on & Hartmami' <199&): (JHK) 
[Skrutskie et al. (2006); (Spitzer IRAQ Robitaille et al. (2007); 
(IRAS )iWeaver & Jones ( 1992); (submm) lAndrews & Wilhamsl 
(l2005h : (mm)£ietuetaL (.2006) : 



With the size of the shadow cast by the inner disk con- 
strained, we proceed by fitting the outer disk. We find a disk 
mass of 0.3 Mq using p=l. This disk mass is much higher 
than p reviously reported, for instance by Andrews & Williams 
(l2005h . This can be attributed to the millimeter opacity of our 
big particles, which is a factor of five lower at 1 .3 mm. Corrected 
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Fig. 3. Radial t= 1 surfaces for our best-fit disk models of LkCa 
15. The solid line shows the disk model with an optically thick 
inner disk. The inner disk casts a shadow on the outer disk, 
and only the disk surface is directly illuminated. The dotted line 
shows the disk model with an optically thin dust shell. The inner 
rim of the outer disk is fully illuminated, and the outer disk is 
much flatter In this plot, photons from the central star travel on 
horizontal lines (constant z/r). 



fo r this, the disk mass fo r LkCa 15 does not differ with from that 
of lAndrews & Wilhamsl (l2005h . 

We lower the outer disk surface by using a sedimentation 
parameter of *I'aii=0.35 and a big particle fraction of 94%. We 
also settle the big grains to the midplane (*I'big=0.5) to fit the 
relatively low flux at 350 fim. A model without stratification 
(*Pbig=l) overpredicts this flux a little, but is still consistent with 
the large errors at that wavelength. The inner rim of the outer 
disk is almost completely shadowed by the inner disk, and only 
the disk surface is illuminated (Fig.[3]l. The reduction in surface 
height compared to a well-mixed model without big part icles is 
four, comparable to that of other T-Tauri stars (.D'Alessio et al.l 
l200a) . 



3.2. Optically thin inner disk 

To investigate the effects of the optical depth of the inner disk 
on the outer disk structure, we first construct a disk model with- 
out an inner disk. In this model, the inner rim of the outer disk 
is fully illuminated, and its contribution to the SED from 25yum 
to lOOyLfm will be much larger compared to a shadowed outer 
disk. To fit the far infrared flux, we need stronger sedimentation 
(*I'aii=0.25) and more big particles (fmn, =99%), so the outer disk 
becomes twice as flat (Fig. [3] dotted line). We find a total reduc- 
tion in the surface height of a factor of eight, significantly higher 
than in other T-Tauri stars, which was also reported by Esp07. 

To fit the complete SED, we follow Esp07 and add an inner 
disk between 0. 1 and 5 AU. If we assume this disk is in vertical 
hydrostatic equilibrium, we fit a dust mass in small grains of 
2 X 10"'" Mq as constrained by the flux in the near and mid 
infrared. This inner disk however is not optically thin to starlight. 
Its radial t=1 surface reaches the same height (z/r~0.1) as the 
outer disk, casting a shadow over the entire outer disk which 
results in a far infrared flux drop by a factor of ten (dotted line 
in the lower panel of Fig. |2]i. 

One way to minimize the effect of this shadow is to make the 
inner disk extremely flat (z/r ~ 0.01) (dashed line in the lower 
panel of Fig. |2]i. Although this model fits the far infrared flux 
reasonably well, the inner disk does not reprocess enough of the 
stellar radiation to reproduce the near infrared excess. Another 
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Model parameters 



Parameter 


Value 


Range 


T TKl 

A eiT L^^J 


4350' 




T FT 1 


0.74' 




M TM^ 1 

ivi* LiviQ J 


0.97' 




d[pc] 


126 


110-140 


Rin,outer [AU] 


46 




Rout,outer [AU] 


800 




Mdisk,outer [Mo ] 


0.3 


0.01-1.0 


P 


1.0 




Rin,iniier [AU] 


0.1 




Optically thick Optically thin 
inner disk dust shell 


l^out,inner [AU] 


1 5 


0.2-5 


Mdust,mner [Mq ] 


4 X 10-* 1 X 10-" 


lo-'^t-io-^ 




1.00 


0.01 - 1.00 




0.35 0.25 


0.15-1.0 


fmm [%] 


94 99 


0-99.99 




0.5 1.0 


0.01-1.0 


A, [mag] 


1.2 0.95 


0.0-2 



Table 1. Stellar and best-fit parameters for the tw o diff erent 
geometries. Refere nces - ( 1) .Kenvon & Hartmann (Il996l) : (2) 
ISimon et al.1 ( 120001) 



way to avoid the shadow is to increase the disk scale height much 
beyond its hydrostatic equilibrium value. By spreading the same 
amount of dust over a larger solid angle, its radial optical depth 
decreases. For the material to become optically thin, the scale 
height has to be so high that we will model it as a spherical 
shell around the star (solid line in the lower panel of Fig.|2|i, but 
the exact geometry is hard to constrain as long as the material 
remains optically thin. We need only 1 x 10"' ' Mq of dust in this 
shell because it does not shield itself and is therefore hotter. It is 
also in our line of sight and reddens the starlight by an additional 
A„=0.25. 



4. Discussion and conclusions 

We have studied the radiative transfer effects in the gapped disk 
of LkCa 15 using MCMax, a self-consistent 2D radiative transfer 
code. We show that only an optically thin dust shell or an opti- 
cally thick inner disk can explain the SED in the near and mid 
infrared. The outer disk geometry differs significantly depending 
on the inner disk properties. With an optically thick inner disk, 
the outer disk must rise out of the shadow cast by the inner disk. 
With an optically thin dust shell, there is no shadow and the outer 
disk becomes twice as flat. 

These models are similar to but not the same as those pro- 
posed in Esp07. The latter do not take into account the shadow 
cast by an optically thick inner disk on the outer disk. They find 
the same z/r for the inner disk rim (0.01 AU / 0. 12 AU =0.08) as 
for the outer disk rim (4 AU / 46 AU =0.08). This means that the 
outer disk wall cannot be illuminated, as it lies in the shadow of 
the inner disk (see Fig. |2] dotted line). The model with an opti- 
cally thin inner disk can be reproduced, with the exception that 
to become optically thin, its scale height needs to be so high that 
it effectively becomes a dust shell. 

Esp08 find the near infrared excess to be a blackbody and 
conclude that the inner disk has to be optically thick. On the 
other hand, Vinkovic et al. (2006) concluded that compact dust 
halos can explain the near infrared bumps in Herbig Ae/Be stars 
and that their existence agrees well with interferometric data. 




Fig. 4. 20 micron images of our best-fit disk models of LkCa 15. 
The left panel shows a model with an optically thick inner disk 
of 1 AU that casts a shadow on the rim of the outer disk. The 
right panel shows a 5 AU optically thin dust shell and a fully 
illuminated outer rim. Note that if the rim of the outer disk is 
shadowed, much more radiation emerges from the disk surface. 

An open question remains as to why the inner disk would 
have so high a scale height that it becomes an optically thin dust 
shell. Such a high scale height for small dust grains could be due 
to a halo of unknown origin, maybe residual infall, but such an 
infall would not lead to a dust cloud confined to <5AU. Another 
possibility could be violent scattering and disruption of planetes- 
imals that could create a cloud of dust with a scale height deter- 
mined by dynamic processes rather than the hydrostatic equilib- 
rium of gas. However, there is currently no additional evidence 
supporting this scenario for LkCa 15. Furthermore an optically 
thick inner disk leads to a model with outer disk settling require- 
ments similar to what is seen in other T Tauri stars. Therefore 
we select the model with the optically thick inner disk as the 
preferred one. 

To distinguish between geometries with a dust disk or a dust 
shell - or anything inbetween - a near-infrared interferometer can 
be used to directly probe the inner regions. But because a shadow 
from the inner disk connects the inner and outer disk geometry, 
imaging of the outer disk can also reveal the nature of the inner 
disk (Fig. |4|. If the inner disk is optically thick, radiation from 
the outer disk will be more extended because the contribution 
from the shadowed inner rim of the outer disk is small, while 
the disk surface has a larger flaring angle and contributes to the 
20 jum flux as well. For an optically thin inner dust shell, the 
emission will be more centrally peaked because the outer disk is 
flatter and its rim fully illuminated. 

We can draw one general conclusion about transitional ob- 
jects: the radial optical depth - and therefore the vertical struc- 
ture - of dust within an inner hole determines the illumination of 
the outer disk and ergo its geometry. This is especially impor- 
tant in the light of observations with the Spitzer space telescope, 
where large nu mbers of disks wit h dust-filled inner holes have 
been identified (iNaiita et alj|2007h . 
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